A new model for the double well potential 
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A new model for the double well potential is presented in the paper. In the new potential, the 
exchanging rate could be easily calculated by the perturbation method in supersymmetric quantum 
mechanics. It gives good results whether the barrier is high or sallow. The new model have many 
merits and may be used in the double well problem. 
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Introduction 

The quantum tunneling in double well potential have 
been extensively applied in many branches of physics. In 
spite of being almost as old as quantum mechanics, the 
tunneling dynamics of the double well are one of the most 
active research areas. For example, they appeared in the 
mean-fields dynamics of Bose-Einstein condensates, the 
recent development of ion trap technology, the ultracold 
trapped atoms theory and applications, etc[l]-[5]. In all 
these the single particle energy splitting, the ener- 

gies and the wave functions in the double well are needed 
for the relevance problem, the usual potential model of 
the double well for the single particle is 



V(x) 



-ax 



bx\ 



(1) 



The model gives rise to the challenge to the mathemati- 
cal and computational techniques. The focus of the prob- 
lem of finding the energy splitting due to the tunneling 
turned out to be difficult. It can not be treated by the 
regular perturbation methods. Then many methods ap- 
peared to attack the problem, like the instanton method, 
the WKB approximation method, Variational method, 
numerical method, and the improved WKB method, the 
super-symmetry quantum method, etc. The validity of 
these methods are still limited. 

It is the form of the potential that make all the ef- 
forts of the above mentioned methods not satisfactory. 
Because the property of the potential (1) at the infinity 
x = oo, the singular point of the Schrdedinger's, changes 
steeply from b = to b ^ 0, it presents one of the great- 
est challenge to the methods of solving differential equa- 
tions. Actually, no matter how small the parameter b 
might be, the perturbation term bx 4 would absolutely 
exceeds the un-perturbation term ax 2 when x is large 
enough. So the property at the infinity of the potential 
and the Schrodingcr equation is mainly determined by 
the term bx , which is no longer the small perturbation 
term. Hence the regular perturbation method for small 
b to solve the double problem does not work for sure. It 



is natural that there still have some controversial conclu- 
sions on how reliable the above mentioned methods may 
be. 

As stated before, it is the potential form of (1) that 
gives rise to the difficulties shown in treating it. The 
best way to the problem is to have some models other 
than that of the (1). The ideal model should be solved 
analytically. There is a model of the delta split harmonic 
oscillator almost meets the requirement [6]: 



1 



V(x) = -muj' i x' i + kS(x) 



(2) 



Though there are some possible situations for the ideal 
or perfect point-like potential, it is somehow still a math- 
ematical idealization. 

Hence another solution that naturally comes forth is 
to give a new physical model. The new model should 
have the traits of the double potential. It better had 
completely solvable property, if not, it should allow the 
regular perturbation methods to apply. In the paper, we 
propose a physical model meeting all these requirements 
for the double well potential. Actually, we have found it 
in studying the spheroidal equations. 

In solving spheroidal wave functions, we have con- 
nected them with the following equation 
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and the boundary conditions for ^ become 

#U=o = #U=, = (4) 
The the potential V(x, a, m) is 
1 
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V(x, a, m) = — — — a cos 2 x H =— 4- (5) 
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where the parameters m are integers. The eqn.(3) have 
been solved by the perturbation methods in supersym- 
metry quantum mechanics [9], [10]. See appendix for the 
detail of the complete solution. 

The potential is symmetry with the center at x = \. 
It has many interesting properties. When the parameter 
m = 0, the potential V(x, a, 0) behaves as the figures 
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FIG. 1: comparison of the potentials V(x,a,0) in the case 
m — 0. The blue, red and yellow lines are the potential with 
a = 0, —0.23, —0.4. The figure shows the potential changes 
its properties when a < — ~ . 



FIG. 2: comparison of the potentials V(x,a,2) in the case 
m = 2. The blue, red and yellow lines are the potential with 
a — 0, 0.6, 5. The figure shows the formation of the double 
well potential when a > 4 — ~ . 



shown for the ranges a > — ^ and a < — j respectively. 
Only in the case a < — j, does it possess two local mini- 
mums and look like a reverse double well potential. 

When m ^ 0, the potential changes greatly, and ap- 
proaches positive infinity at x = 0, tt. In the range 
a < m 2 — |, the difference of the potentials between 
a = and a ^ is artificial. The fundamental changes 
come when m 2 — \ < a. The nice property of the po- 



tential comes for 



< a: it resembles that of the 



form of the double potential. The potential has two equal 
minimum values 



Vrr, 



a(m 2 -\)~\a 



at the symmetrical points x m i nt i, x m i n ^ satisfying 
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The maximal value is achieved at the center point x 
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FIG. 3: comparison of the potentials V(x, a, 0) in the case 
m = 2. The blue, red lines are the potential with a = 10, 5. 
The figure shows the double well depth decreases as a de- 
creases. 



the barrier's height is 



V, 



max v 771 tri 
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(9) 



It is easy to see that the depth of the double well is com- 
pletely determined by the parameter a when the integral 
quantity m is fixed. The depth will increase as a goes 
up. Small a corresponds the sallow double well. 

The great advantage of the potential (5) lies in that the 
properties of the singularities of the potential do not be 
influenced by the term —a cos 2 a;, which only adds a con- 
stant —a to the infinity at x — or x = tt. This stands 



in contrast with that in the potential (1). The term 
—a cos 2 a; actually brings the double- well-form potential 
into reality. This is just what we want to achieve. Hence, 
one can treat the problem by the usual perturbation 
methods, which have been studied in series papers [7]- 
[11]. 

In the double well problem, the tunneling rate is the 
energy difference t = E\ — Eq between the lowest wave 
functions. The focus is how accurate one can get the 
tunneling rate through the double barrier. In the new 
model, this corresponds to the cigen-encrgy difference of 
t = E\. m — E(j. m , which has been given in previous works 



[9], [10] 



Ei vm - Eo-,m = R(ai), i?(ai) = ^V-R^oi), (10) 



a\ = (Aij), Ai^ = 1, 

«2 = \Bi,j), Bi_j = Bi ! j(Ai i k) 

R Q = (2m + 2) 

B hl +A hl 
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(11) 
(12) 
(13) 

2m + 3 (W) 
R n {A. u ) = [c+^Anj] - c-^Bnj)] ,n>2 (15) 



where 



<i(A*,i) = - [(2m + l)A ,o - 1] A ra>1 a„,i (16) 
C~l(Bn,j) = - [(2m + l)B ,o + 1] Bn.lOn.l (17) 

See the reference [10] for the determination of the quan- 
tities Ojj, Bij. Here we give some example, when m = 2, 



i?o = 12, Ri = -0.19047619, 
R 2 = -0.001374287. 

From the relation (10), we have 



(18) 
(19) 



a = 5->t = E 1 - E Q = 11.01326187, (20) 
a = 10 -> t= 9.9578094 (21) 

respectively. The potential (5) represents the symmetry 
double potential, and the asymmetry potential could be 
represented by 



V(x, (3, m, s) 



(m + scosx) 2 - \ 2 2 
— — p cos X 
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(22) 



This potential appears in the study of the spin-weighted 
spheroidal functions [15], [16], [17]. The spin-weighted 
spheroidal functions come from the study of the pertur- 
bation problem of the Kerr black hole[12]-[14]. The sim- 
ilar results about the cigcn-valucs and cigcn-functions to 
that of the references [9], [10] will be reported in the fu- 
ture. This potential approaches infinity at x = 0, tt and 
might have two minimums and one maximum in some 
range of the parameter a. The graph of the potential 
is similar to that of the asymmetry double well in some 
way: their distances to the center x = ^ not equal. How- 
ever, the two depthes of the asymmetry potential are not 
equal. 



The double well potential (1) have been applied ex- 
tensively in physics. Their fundamental properties are in 
abundance depending the parameters a, b. In the pa- 
per, a new type of potential (5) is proposed. They are 
shown to have the similar fundamental properties as that 
for the double well ones (1). The nice properties of the 
new potential is that their Schrodingcr equations could be 
solved by the usual perturbations methods, which stands 
in vivid contrast against the one (1). Furthermore, their 
solutions in perturbation series of the parameter a are 
already given in the [9], [10]. All these results are useful 
for their application in the double well problem. 
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